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Abstract. Breather solutions of the modified Korteweg-de Vries equation are shown to be 
globally stable in a natural H 2 topology. Our proof introduces a new Lyapunov functional, at 
the H 2 level, which allows to describe the dynamics of small perturbations, including oscilla- 
tions induced by the periodicity of the solution, as well as a direct control of the corresponding 
instability modes. In particular, degenerate directions are controlled using low-regularity con- 
^s^J servation laws. 

I—* 

1. Introduction 

This paper deals with the nonlinear stability of breathers of the focusing, modified Korteweg-de 
^ Vries (mKdV) equation 

u t + (u xx + u 3 ) x = 0. (1.1) 

Here u = u(t,x) is a real-valued function, and (t, x) G M 2 . The equation above is a well known 
completely integrable model [131 HI [2D], with infinitely many conserved quantities, and a suitable 
Lax-pair formulation. The Inverse Scattering Theory has been applied by many authors in order 
to describe the behavior of solutions in generality, see e.g. [TJ |20j and references therein. 



(N 



Solutions u(t,x) of (1.1) are invariant under space and time translations, and under suitable 



scaling properties. Indeed, for any to,xo € K, and c > 0, both u(t—to, x—xq) and c 1 / 2 w(c 3 / 2 t, c x l 2 x) 



are solutions of ( 1.1). Finally, if u(t, x) is a solution of (1.1), then u(—t, —x) and — u(t, x) are also 
solutions. 



> 

in 

CO 

VO M[u] (t) := \ [ u 2 (t,x)dx = M[«](0), (1.2) 



On the other hand, standard conservation laws for (1.1) at the .H^-level are the mass 



and energy 



< £?[«](*) :=i / u 2 x (t,x)dx-~ f u A (t,x)dx = E[u}(0). (1.3) 

A satisfactory Cauchy theory is also present at such a level of regularity or even lower, see e.g. 
Kenig-Ponce-Vega [TS], and Colliander et al. [11". From the Inverse Scattering Theory, the 
evolution of a rapidly decaying initial data can be described by purely algebraic methods. Solutions 
are shown to decompose into a very particular set of solutions (see Schuur |33j). described in detail 
below. 



Indeed, equation (1.1) is also important because of the existence of solitary wave solutions 



called solitons. These profiles are often regarded as minimizers of a constrained functional in the 



H 1 -topology. For example, mKdV (1.1) has solitons of the form 

u(t,x) = Q c (x - ct), Q c (s) := y/cQ(y/cs), c> 0, (1.4) 

with 

/2 

Q(s) := \ . . = 2\/29Jarctan(e s )]. 
cosh(s) 
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By replacing (1.4) in (1.1), one has that Q c > satisfies the nonlinear ODE 

Q"-cQ c + Q 3 c =0, Q c eH\M). 



(1.5) 



Moreover, as a consequence of the integrability property, these nonlinear modes interact elastically 
during the dynamics, and no dispersive effects are present at infinity. In particular, even more 
complex solutions are present, such as multi-solitons (explicit solutions describing the interaction 



of several solitons [H]). For example, the 2-soliton solution of ( 1.1 ) is given by the four-parameter 
family 

e Sl + e S2 



U 2 ■= U 2 (t,x;c 1 ,C2,x 1 ,x 2 ) = 2\fld x 



arctan 



1 



with si 



ci(x—cit) +xi, s 2 := y/c&(x - c 2 t) +x 2 , and p := Here ci, c 2 > 0, ci ^ c 2 , 

are the corresponding shift parameters. In particular, 



are the associated scalings, and xi,x 2 G 
U 2 satisfies 

lim \\U 2 (t) - Q C1 (- - cit - xf) 

t— >±oo 

for some given x^ g K, depending on (c\, c 2 ). 



Qc 2 (- ~ C 2 t- X 2 )||ifi( R ) = 0, 



The study of perturbations of solitons and multi-solitons of (1.1) and more general equations 



leads to the introduction of the concepts of orbital, and asymptotic stability. In particular, since 
energy and mass are conserved quantities, it is natural to expect that solitons are stable in a 
suitable energy space. Indeed, H 1 -stability of mKdV and more general solitons and multi-solitons 
has been considered e.g. in Benjamin [8 , Bona-Souganidis-Strauss jTUj, Weinstein [37] . Maddocks- 
Sachs [12], Martel-Merle-Tsai [27] and Martel-Merle [23 HI]- instability of KdV solitons and 
multi-solitons has been proved in Merle- Vega [25] and Alejo-Muhoz-Vega [BJ. On the other hand, 
asymptotic stability properties have been studied by Pego- Weinstein [21] and Martel-Merle [231 

EHESJ- 

One of the main ingredients of the stability argument employed in some of the previous works 
is the introduction of a suitable Lyapunov functional, invariant or almost invariant in time and 
such that the soliton is a corresponding extremal point. For the mKdV case, this functional is 
given by 

H[u](t) = E[u]{t) + cM[u]{t), (1.6) 



where c > is the scaling of the solitary wave, and E[u], M[u] are given in (1.2|-(1.3). A simple 
computation shows that for any z(t) 6 H 1 (M.) small, 



H[Q C + z]{t) = H[Q C ] + f z(Q" c - cQ c + Q 3 C ) 



Q(t) + 0(\\z(t)\\ 3 Hl 



(1.7) 



The first term above is independent of time, while the second one is zero from (1.5). It turns out 



that the third term Q(t) is positive definite modulo two directions, related to the invariance of the 
equation under shift and scaling transformations (see the second paragraph above). Modulation 
parameters are then introduced in order to remove those instability modes. Once these directions 



are controlled, the stability property follows from (1.7) 



In addition to the special solutions mentioned above, there exists another nonlinear mode, of 
oscillatory character, known in the physical and mathematical literature as the breather solution, 
and which is a periodic in time, spatially localized real-valued function. Indeed, the following 
definition is standard (see [351 HO] and references therein): 



Definition 1.1. Let a, 



B a ,p{t, x) :=2v / 2(?£c arctan 



The breather solution of mKdV (1.1 1 is explicitly given by 



/ P sin(a(a; + 8t)) \ 
\a cosh(/3(a; + ^t)) J 



=2^2/3 sech(/3(x + 7i)) 



a cosh(/3(.T + 7^)) 

cos(a(x + St)) - (0/a) sin(a(a; + 6t)) tanh(/3(x + 7*)) 



1 + (/3/a) 2 sm 2 (a{x + 5t)) sech^/^rr + jt)) 



(1.8) 
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with 

2 oo2 o„2 



8:= a* -3/3* 7 := 3c^ - ^. (1.9) 

Note that breathers are periodic in time, but not in space, and this will be essential in our 
proof. A simple but very important remark is that 5 ^ 7, for all values of a and /3 different from 
zero. This means that variables x + St and x + ft are always independent. Indeed, if 6 = 7, one 



has from (1.9 1 2(a 2 + f3 2 ) — 0, which means a = /3 = 0, a contradiction. 

Additionally, note that for each fixed time, the mKdV breather is a function in the Schwartz 
class, exponentially decreasing in space, with zero mean: 

B a ,ji = 0. 

Moreover, from the scaling invariance, one has c 1 / 2 B a) p{c 3 / 2 t, c x l 2 x) = B c i/2 a c i/2o(i, x), for all 
c > 0, and B_ a p = B a p, B a _p = —B a p. Therefore, we can assume a,/3 > 0, with no loss of 
generality. Finally, we will denote j3 and a as the first and second scaling parameters, and —7 
will be for us the velocity of the breather solution. 



For the sake of completeness, we briefly comment the two limits /3/a -C 1 and a = in (1.8). 
The first one allows to simplify the expression for the breather to 

B a ,p(t,x) w 2v / 2/3cos(a(x + 5t))sech(/3(x + ~/t)) + o(-), 

and from a qualitative point of view, it shows explicitly its wave packet nature, as an oscillation 
modulated by an exponentially decaying function (see e.g. Fig. [I]). The second case is obtained 
by formally taking the limit a — > in ( |1.8| , 

B oA t,x) -^[arctan(^-^)]. (1.10) 

This is the well known double pole solution of mKdV (see e.g. [3D]), which represents a soliton- 
antisoliton pair traveling in the same direction and splitting up at logarithmic rate. 




Figure 1. Left: mKdV breather (1.8) with a = 9,0 = 1 at t = 0. Right: mKdV 



breather (1.8) with a = 1.5,0= 1 at t = 0. 



Note that from the invariance under space and time translations, given any to,xo € M, the 
function B a p(t — to,x — Xq) is also a breather solution. This fact allows to define a four-parameter 
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family of solutions 

B a ^(t,x;xi,x 2 ) := B a ^(t -t ,x- x ) = 2\f2d x 



arctan 



fP sin(gyi) \ 
Va coshf/Juo) / 



with yi := x + St + Xi, y 2 :— x + ft + x 2 , 

X\ - x 2 



to : = 



2(a 2 + /3 2 ) 

Note that from this formula one has, for any k E Z, 

kn 



and xq := 



■ a cosh(^y 2 )- 
8x2 - fX\ 



2(a 2 +/3 2 



B a »(t,x;xi H ,a; 2 ) = (-1) 'B a ^(t,x;xi,x 2 ), 

a 



(1.11) 
(1.12) 
(1.13) 



which are also solutions of (1.1). These identities reveal the periodic character of the first trans- 
lation parameter. 



In the same way, from (1.10) one can define a three-parameter family of double pole solutions 
Bo,p{t; xi,x 2 ), with xi,x 2 e R. 



a = 7 
0=2 




Figure 2. Evolution of the mKdV breather (1.8) with a = 7,/3 = 1 at instants 
£ = —0.09, t = 0, and t — 0.13. Note that with the selected values of a,f3, the 
velocity is given by 7 = 3a 2 — ft 2 = 146 > and then the breather moves to the 
left. 



Let us come back to breather solutions. We claim that they can be formally associated to 
the well known mKdV 2-solitons. Indeed, they have a four-parameter family of symmetries: 
two scaling and two translations invariances (note that the equation that we consider is just one 
dimensional in space). However, unlike 2-soliton solutions, breathers have to be considered as fully 
bounded states, since they do not decouple into simple solitons as time evolves. Another intriguing 
fact is that, as far as we know, breathers are only present in some very particular integrable 
models, such as mKdV, NLS and sine Gordon equations, among others. 

Let us recall now some relevant physical and mathematical literature. From the physical point of 
view, breather solutions are relevant to localization-type phenomena in optics, condensed matter 
physics and biophysics [7 ]. In a geometrical setting, breathers also appear in the evolution of 
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closed planar curves playing the role of smooth localized deformations traveling along the closed 
curve [2]. Moreover, it is interesting to stress that breather solutions have also been considered 
by Kenig, Ponce and Vega in their proof of the non-uniform continuity of the niKdV flow in 
the Sobolev spaces H s , s < \ [19] , On the other hand, they should be essential to completely 
understand the associated soliton-resolution conjecture for the mKdV equation, according to the 
analysis developed by Schuur in (33] . An essential problem in that direction is to show whether 
or not breather solutions may appear from general initial data, and for this reason to study their 
stability is the fundamental question. Numerical computations (see Gorria-Alejo-Vega [3]) show 
that breathers are numerically stable. However, the simple question of a rigorous proof of orbital 
stability has become a long standing open problem. 

In this paper, we give a first, positive answer to the question of breathers stability. Our main 
result is the following 

Theorem 1.2. mKdV breathers are orbitally stable in their natural H 2 -topology. 



A more detailed version of this result is given in Theorem 6.1 As we will see from the proofs, 
the space H 2 is required by a regularity argument and by the very important fact that breathers 
are bound states, which means that there is no mass decoupling as time evolves. However, our 
argument is general and can be applied to several equations with breather solutions, and moreover, 
it introduces several new ideas in order to attack the stability problem in the energy space. In 
addition, our proof corroborates, at the rigorous level, some deep connections between breathers 
and the 2-solitons of mKdV. 

Let us explain the main steps of the proof. First, we prove that breathers satisfy a fourth-order, 



nonlinear ODE (equation (3.6)). The proof of this identity is involved, and requires the explicit 
form of the breather, and several new identities related to the soliton structure of the breather. It 
seems that this equation cannot be obtained from the original arguments by Lax [21j . since the 
dynamics do not decouple in time. Our second and more important ingredient is the introduction 
of a new Lyapunov functional (see (5.2 1), well-defined in the H 2 topology, and for which breathers 



are surprisingly not only extremal points, but also local minimizers, up to symmetries. This 
functional also allows to control the perturbative terms and the instability directions that appear 



during of the dynamics, the latter as consequences of the symmetries described by ( 1.8 1. From the 
proofs, we will see that breathers have essentially three directions of instability, two associated to 
translation invariances, and a third one consequence of the particular first scaling parameter /3. In 
order to prove that there is just one negative eigenvalue, we make use of a direct generalization of 
the theory developed by L. Greenberg [TJ], which deals with fourth order eigenvalue problems. We 
then modulate in time in order to remove the spatial instabilities. This is an absolutely necessary 
condition in order to obtain an orbital stability property. However, we do not modulate the scaling 
instabilities. Instead, we control the dynamics first replacing the corresponding negative mode by 
a more tractable direction, the breather itself, and using the mass conservation law. This technique 
was first introduced by Weinstein in [36] . A very surprising fact is that the so-called second scaling 
parameter, associated to oscillations, is actually a positive direction when enough regularity is on 
hand, and even if it has an L 2 -critical character. 

Our functional is reminiscent of that appearing in the foundational paper by Lax |21j , concerning 
the 2-soliton solution of the KdV equation, 

u t + (u xx + u 2 ) x = 0, 

and generalized to the KdV TV-soliton states by Maddock-Sachs [22] . This idea has been success- 
fully applied to several 2-soliton problems, for which the dynamics decouples into well-separated 
solitons as time evolves, see e.g. Holmer-Perclman-Zworski [IB], Kapitula [T7j, and Lopes-Neves 
[29] . for the Benjamin-Ono equation. However, there was no evidence that this technique could 
be generalized to the case of even more complex solutions, such as breathers. Compared with 
those results, our proofs are more involved, and computations are sometimes a nightmare. We 
have preferred to split the proof of the main theorem into several simple steps. 
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We believe that our result can be improved to reach the H 1 level of regularity, but with a 
harder proof. It seems clear that a better understanding of the H 1 dynamics requires a detailed 
study of modulations on the scaling parameters. In particular, the Martel-Merle-Tsai technique 
[27] seems to fail in this case due to the absence of a clearly decoupled mass dynamics. One can 
also consider a suitable asymptotic stability property, in the spirit of |23j . However, note that 
the Martel-Merle [231 1211 results are difficult to generalize to the current case of study since 
breathers may have negative velocity, and therefore they can interact with the linear part of the 
dynamics. We conjecture that breathers are asymptotically stable in the case of positive velocities. 



Remark 1.1. The methods employed in the proof of Theorem 1 1 . 2 1 seem do not apply in the limit 
a — > 0, which is expected to be unstable, according to the numerical computations performed by 
Gorria-Alejo-Vega g]. 

Remark 1.2. The natural complement of our study is to consider the sine Gordon equation 

Utt — u xx + sin u — 0, u(t, x) € R. 

Since this integrable equation has also breather solutions (see e.g. Lamb |20) ) . we expect similar 
results, but with more involved proofs at the level of the linearized problem (we deal with matrix 
operators). Indeed, following the present proof, we can guess that sine-Gordon breathers are 



H 2 x H 1 stable provided Lemma 4.3 and Proposition 4.8 hold for the associated spectral elements. 



Additionally, the focusing Gardner equation 

u t + (u xx + u 2 + [iu 3 ) x = 0, n > 0, 



is the natural generalization of (1.1). In particular, it has a family of breathers indexed by the 
additional parameter fi (see [3"2"ll3"]V We expect to consider some of these problems in a forthcoming 
publication (see [5])- 

In a more qualitative aspect, we think that our results are in some sense a surprise, because 
any nontrivial perturbation of an integrable equation with breathers solutions should destroy the 
existence property. Several results in that direction can be found e.g. in GH] and references 
therein (for the sine Gordon case). Those results and the present paper suggest that stability is 
deeply related to the integrability of the equation, unlike the standard gKdV TV-soliton solution 
[27]. 

Finally, let us explain the organization of this paper. In Section[2]we study generalized Weinstein 
conditions satisfied by breather solutions. In Section [3] we prove that any breather profile satisfies 
a fourth order, nonlinear ODE. Section [4] is devoted to the study of a linear operator associated 
to the breather solution. In Section [5] we introduce new Lyapunov functional which controls the 
dynamics. Finally, in Section[6]we prove a detailed version of Theorem 1.2 



Acknowdlegments. We would like to thank Yvan Martel, Frank Merle, Carlos Kenig and Luis 
Vega for many useful comments on a first version of this paper. 



2. Stability tests 

The purpose of this section is to obtain generalized Weinstein conditions for any breather B. 
Indeed, for the case of the mKdV soliton (1.4), the mass (1.2) and the energy (1.3| are given by 
the quantities 



M[Q C ] = ~ 



Q 



2 - 2c 1 ' 2 



E[Q C 



r?l 2 E[Q\ = -\<?'^M[Q\ = -\c*l 2 



< 0. 



(2.1) 



(2.2) 



3 L " J 3" 

These two identities show the explicit dependence of the mass and the energy on the soliton scaling 
parameter. In particular, the Weinstein condition [37] reads, for c > 0, 

d c M[Q c ] = c- 1/2 > 0. (2.3) 
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This condition ensures the nonlinear stability of the soliton. We consider now the case of mKdV 
breathers. Surprisingly enough, the mass of a mKdV breather only depends on the first scaling 
parameter (3. In other words, it is independent of a. 

Lemma 2.1. Let B = B a ^ be any mKdV breather, for a, j3 > 0. Then 

M[B){t) = 2pM[Q] = 4/3. (2.4) 

Proof. We start by writing the breather solution in a more tractable way. From the conservation 
of mass and invariance under spatial and time translations, we can assume X\ = x 2 = t = in 
( |1.11[ ). We have then 

'cos(ax) — {13 /a) sin(ax) tanh(J3x) i 2 



B 2 (0,x) = 8/3 2 sech 2 (/3x) 



1 + (/3/a) 2 sm 2 (ax) sech 2 (/3x) 
Expanding the square in the numerator, we get after some simplifications 
B 2 (0,x) = 8a 2 (3 2 x 

'a 2 cosh 2 (/3a;) cos 2 (ax) + /3 2 sin 2 (ax) sinh 2 (/3x) — 2a/3sin(ax) cos(ax) sinh(/3x) cosh(/3x) 



(a 2 cosh 2 (/3x) + f3 2 sin 2 (ax)) 2 

Now the purpose is to use double angle formulas to avoid the squares. More precisely, it is well 
known that 

cos 2 (ax) = -(1 + cos(2aa;)), sin 2 (aa;) = - (1 — cos(2ax)), (2.5) 



and 



cosh 2 (/3a;) = ^(1 + cosh(2/3ir)), sinh 2 (^) = i(cosh(2/3a;) - 1). (2.6) 



We replace these identities in the previous expression above. We obtain 

8a 2 /3 2 h a ^(x) 



B 2 (0,x) 



(a 2 + (3 2 + a 2 cosh(2/3a;) - (3 2 cos(2ax)) 2 ' 
with 

h a .j)(x) := a 2 - /3 2 + (a 2 + /3 2 )(cos(2ax) + cosh(2/3a;)) 

+ (a 2 - /3 2 ) cos(2ax) cosh(2/3x) - 2a/3 sin(2ax) sinh(2/3a;). (2.7) 

In what follows, let 

fa,p(x) : = a 2 + /3 2 + a/3 sin(2aa;) - f3 2 cos(2ax) + a 2 (smh(2/3x) + cosh(2/3a;)), 

and 

9a,d(x) := a 2 + (3 2 + a 2 cosh(2/3x) - /3 2 cos(2ax). 

It is clear that 

fa b( x ) = 2af3[a cos(2aa;) + (3 sin(2ax) + a cosh(2/3ir) + a sinh(2/3a;)], 

and 

g' aJ3 (x) = 2a(3[f3 sin(2ax) + asinh(2/3x)]. 
Therefore, after a lengthy but direct computation, 

f a A x )9cA x ) - fa,fs{x)g' a ^(x) = 2a 2 [3h a ^(x), 

and then 

B 2 (0,x) = = W^Y. 

In conclusion, we have proved that 

1 f x „ 2 . , , 2/3[a 2 + ^ 2 + a/3sin(2aa;)-/3 2 cos(2aa;) + a 2 (sinh(2/3a;)+cosh(2/3a;))] 
±> (0, s)as — 



2 7_ 00 v ' ' a 2 + /3 2 + a 2 cosh(2/3ir) - ^ 2 cos(2ax) 

Taking limit as x — > +oo, we get the desired conclusion. □ 
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Remark 2.1. Note that the last integral above does not change if we consider a general breather, 
of the form ( |1.11[ ) . Indeed, our proof does not require the time independence of the solution. Then 
we get 

i r x 



M a! p{t,x) 



2 

2/3 [a 



— CO 

2 



/3 2 + a/3 sin(2ayi) - (3 2 cos(2ayi) + a 2 (sinh(2/3y 2 ) + cosh(2/3y 2 ))] 



a 2 + f3 2 + a 2 cosh(2/3y 2 ) - P 2 cos(2ayi) 
with yi and y 2 defined in ( |1.11[ ). This last expression will be useful in Lemma 2.3 



, (2.8) 



A direct consequence of the results above are the following generalized Weinstein conditions: 

Corollary 2.2. Let B = B a p be any mKdV breather of the form ( |1.11[ ). Given t € K fixed, let 

A a B:=d a B, and ApB := dpB. (2.9) 

Then both functions A a B and ApB are in the Schwartz class for the spatial variable, and satisfy 
the identities 

(2.10) 
(2.11) 



and 



d a M[B] = J BA a B = 0, 

dpM[B] = I BApB = 4 > 0, 
Jm 



independently of time. 



Proof. By simple inspection, one can see that, given t fixed, A a B a .p and ApB a .p are well-defined 
Schwartz functions. The proof of (2.10) and (2.111 is consequence of (2.4|, and the definition of 
mass Ol). □ 



Remark 2.2. Comparing (2.101 and (2.111 with the Weinstein condition (2.3), we may think that 
the second scaling parameter a is L 2 -critical. On the opposite side, the first scaling f3 can be seen 
as a stable parameter. 

Lemma 2.3. Let B = B a jj be any breather of the form ( |1.11[ ) , with a, f3 > 0. Then we have 
(1) B = B x , with B = B a ^ given by the smooth L°° -function 



B(t,x) := 2V2arctan (- S "?^V )- 
\a cosh(/??/ 2 ) / 

(2) For any fixed t £ K, we have B t well-defined in the Schwartz class, satisfiying 



B . 



B t + B 



0. 



(3) Finally, let M a ,/3 be defined by (2.8). Then 

1 



B 



-B 4 + 2BB t - 2{M a 



(2.12) 



(2.13) 



(2.14) 



Proof. The first item above is a direct consequence of the definition of B a ^ in ( |1.11[ ). On the 
other hand, (2.13) is a consequence of (2.12) and integration in space (from — oo to x) of (1.1). 
Finally, to obtain (2.14) we multiply (2.13) by B x and integrate in space. □ 



Remark 2.3. The reader may compare ( 2.13 1-( 2.14 1 with the well known identities for the soliton 
solution of mKdV: 

1 



Qc ~ cQ c + Ql 



0. 



:Qt = 0. 



We compute now the energy of a breather solution. 
Lemma 2.4. Let B = B a ^ be any mKdV breather, for a, j3 > 0. Then 

E[B] = 2/3(3a 2 - /3 2 )\E[Q}\ = 2fa\E[Q]\. (2.15) 
Let us remark that the sign of the energy is dictated by the sign of the velocity 7. 
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Proof. First of all, let us prove the following reduction 

E[B]{t) = I [ {M a ,p) t (t,x)dx. 



Indeed, we multiply (2.13) by B a p and integrate in space: we get 

13; = / BB, + / B\ 



On the other hand, integrating (2.14), 



B 



1 



B 4 + 2 BB t — 2 / (M a ,p) t = 



From these two identities, we get 



and therefore 



B 4 = 



(M a ,p) f -2 BBt 



[ Bl = ~ / (M a ,p) t - f BB t 
Jr. 6 Jr Jr 



Finally, replacing the last two identities in (1.3), we get (2.16), as desired. 



Now we prove (2.15). From (2.8) and similar to the proof of Lemma 2.1 we have 



{M a ,p) t = 2/3- 



a 2 



(2.16) 



(2.17) 



where, with a slight abuse of notation, f a p and g ai p are given now by 

f a /3 = a 2 + f3 2 + a(3sm(2ayi) - (3 2 cos(2aj/i) + a 2 (sinh(2/3j/ 2 ) + cosh(2/3y 2 )), 

and 

g a ,p = a 2 + P 2 + a 2 cosh(2/3y 2 ) - (i 2 cos(2m/i). 

It is clear that 

(fa,p)t — 2a(3[aS cos(2ayi) + (36 sin(2ayi) + 07 cosh(2/3y 2 ) + a7sinh(2/3j/ 2 )], 

and 

(9a,p)t = 2a/3[/?<5sin(2cn/i) + aj sinh(2/3y 2 )]. 

Therefore 

(fa,p)t9a,P - fa,p{ga,p)t = 2ct 2 f3h a ^, 

and 

h ai p := -fa 2 - 5f3 2 + (a 2 + (3 2 )(5 cos(2a yi ) + j cosh(2f3y 2 )) 

+ {8a 2 - 7/? 2 ) cos(2ayi) cosh(2/3y 2 ) - af3(S + 7) sin(2ayi) sinh(2/3y 2 ). 

In conclusion, 

3 ( a ' 0)t ~ 3(a 2 + f3 2 + a 2 cosh(2£y 2 ) - (3 2 cos(2ayi)) 2 ' 

Now we split h a p into two pieces, according to the parameter 7. From the definition of 7, S and 
h a ,p, we have 

h a ,p = ~f h a,p + 2(a 2 + j3 2 )h a ^ 7 

where 

h ay p := 1 - (a 2 + (3 2 ) cos(2m/i) - a 2 cos(2ay 1 ) cosh(2/3y 2 ) + a/3sin(2m/i) sinh(2^y 2 ). 



(2.18) 



(2.19) 



Note that from Lemma 2.1 more precisely (2.8), 

4a 2 /3 2 jh a ,p 



3(a 2 + f3 2 + a 2 cosh(2/%>) - /3 2 cos(2ayi)) 2 



; 0j = 2^\E[Q]\, 
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then, in order to conclude, from (2.19) and (2.161 we reduce to prove that 

0. 



f 

■J R 



h 



(a 2 + (3 2 + a 2 cosh(2/3y 2 ) - (3 2 cos(2a 2 / 1 )) 2 



(2.20) 



We prove this identity, noting that (^- sin(2aj/i)) x = cos(2aj/i), and from (2.17), 

5a,/3(^sin(2ayi)) 2; - ^- sm(2ay 1 )(g a ^) x = 
la la 

= cos(2ay 1 )[a 2 + [3 2 + a 2 cosh(2/3y 2 ) - P 2 cos(2ayi)] - /?sin(2ayi ) [/3 sin(2ayi) + a sinh(2/3y 2 )] 
= -/3 2 (cos 2 (2ay!) + sin 2 (2ayi)) + (a 2 + (3 2 ) 003(20;?;!) + a 2 003(20;?/!) cosh(2/3y 2 ) 
— a/3 sin(2ayi) sinh(2/3y 2 ) 

= —h a ,p- 
Therefore, 

-sin(2a?/i) 



l.h.s. of (2.20) = lim 



X-S-+OC 2a(a 2 + ft 2 +a 2 cosh(2f3y 2 ) - /3 2 cos(2ayi)) 



0. 



□ 



Remark 2.4. Note that we could follow the approach by Lax .21_, pp. 479-481] to obtain reduced 
expressions for the mass and energy of a breather solution. However, the resulting terms are 
actually harder to manage than our direct approach. 



Corollary 2.5. Let B = B a p be any mKdV breather. Then 

d a E[B] = Uap\E[Q\\ > 0, d p E[B\ = 6(o 2 - (3 2 )\E[Q}\. 



(2.21) 



3. Nonlinear stationary equations 



The objective of this section is to prove that any breather profile satisfies a suitable stationary, 
elliptic equation. 

Lemma 3.1. Let B = B a ^ be any mKdV breather. Then, for all t € K, 

B xt + 2(M a ,p) t B = 2(/3 2 - a 2 )B t + (a 2 + (3 2 ) 2 B. (3.1) 
Proof. We make use of the explicit expression of the breather. An equivalent for the quantity 



has been already computed in Lemma 



and (1.11 1, and using double angle formulas in t 

4V2a(3 



2.4 



see (2.19). On the other hand, from (2.12| 



B, 



le denominator, we have, 
~aS cos(ayi) cosh(/3?y 2 ) - /3~fsm(ayi) sinh(/3jy 2 ) 



From (1.8) we have 



B = 2V2af3 
= A^2af3 



a 2 + (3 2 + a 2 cosh(2/3?/ 2 ) - /3 2 cos(2oyi) 
acos(ayi) cosh(/3?/ 2 ) - j3 sin(ayi) sinh(j8y 2 ) 



a 2 cosh 2 (/3y 2 ) + (3 2 sin 2 (ay 1) 
a cos(ayi) cosh(/3?/ 2 ) — /3 sin(a?yi) sinh(/3?/ 2 ) 
a 2 + fi 2 + a 2 cosh(2/??/ 2 ) - (3 2 cos(2ay x ) 



Therefore, from (3.3) and (2.19), 



(M a ,ff) t B = 



16y/2a 3 P 3 8{t,x) 



{a 2 + f3 2 +a 2 cosh(2(3y 2 ) - P 2 cos(2ayi)) 3 



where, with the definition of h in (2.18), 



d(t,x) := h a ^(a cos(o?/i)cosh(^y 2 ) - fi sin(ayi) sinh(/3y 2 )). 



(3.2) 

(3.3) 
(3.4) 
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Let us compute B xt . First we have from ( |3.3[ ), 

x 



B x = 



(a 2 + /3 2 + a 2 cosh(2/3y 2 ) ~ P 2 cos(2a yi )) 2 ' 



where 



hi(t, x) := -{a 2 + P 2 ) cosh(/3y 2 ) sin(ayi)[a 2 + /3 2 + a 2 cosh(2/3y 2 ) - P 2 cos(2ayi)] 

- 2a/?[acos(ayi) cosh(^j/ 2 ) - P sin(ayi) sinh(/3y 2 )] [/3 sin(2cn/i) + a sinh(2/?y 2 )]. 



Then, 



B x t 



(a 2 + /3 2 + a 2 cosh(2/3 2 / 2 ) - /3 2 cos(2ayi)) 3 1 



where, using g a ^ previously defined in (2.171, 

h 2 (t,x) := -4a/3[^5sin(2ayi) + cry sinh(2/?y 2 )] h 2 \ + g a ,ph 22 , 

and 

h 21 (t, x) := -{a 2 + /3 2 ) cosh(/3y 2 ) sin(aj/ 1 )g a , l g 

— 2a/3(acos(a?/i) cosh(/3y 2 ) — j3 sin(m/i) sinh(/3y 2 ))(/3 sin(2az/i) + a sinh(2/3y 2 )), 

/i 22 (t, a;) := 3 Q:/3 x ( - aS(a 2 + /3 2 ) cos(ayi) cosh(/3y 2 ) - /?7(a 2 + /3 2 ) sin(ayi) sinh(/3y 2 ))g Q ^ 

— 2a^(a 2 + /3 2 )sin(a?; 1 )cosh(/3y 2 )(Q!7sinh(2/3?/ 2 ) + 5f3 sin(2ay 1 )) 

— 2a/3( — a 2 S sin(ayi) cosh(f3y 2 ) + af3-f cos(ayi) sinh(/3y 2 ) 

— f3a5 cos(aj/i) sinh(/3y 2 ) — /3 2 7 sin(a?/i) cosh(/3y 2 )) x (a sinh(2/3y 2 ) + /3 sin(2ayi)) 

— 4a 2 /3 2 (acos(a?/i) cosh(/3y 2 ) — /3sin(ayi) sinh(/3y 2 )) x (7 cosh(2/3y 2 ) + <5cos(2ayi)) 



Then, 



B xt + 2(M a:/s ) t B 



4y/2a(3[h 2 +8a 2 /3 2 



(a 2 +/3 2 +a 2 cosh(2/3y 2 ) - f3 2 cos(2ayi)) 3 ' 
and recalling that 

91,0 = (" 2 +P 2 +a 2 cosh(2/3y 2 ) - ^ 2 cos(2m/i)) 2 

= (a 2 + ,3 2 ) 2 + (a 4 cosh(2/3j/ 2 ) 2 + ^ 4 cos(2ayi) 2 - 2a 2 /3 2 cosh(2/3y 2 ) cos(2ayi)) 
- 2(a 2 + P 2 ){a 2 cosh(2/3y 2 ) - /3 2 cos(2m/i)), 



(3.5) 



collecting terms in (3.5) and taking into account (3.4) and (3.2), after some calculations we have 

rhs of ((351) = 4 ^ x 

' " i - J (a 2 + /3 2 + a 2 cosh(2/3y 2 )-/3 2 cos(2ayi)) 3 

{gl,p x 2(/3 2 - a 2 )(<5acos(ayi)cosh(/3y 2 ) - 7/3sin(ayi) sinh(/3y 2 )) 
+ (a 2 + /3 2 ) 2 (a cos(ayi) cosh(/3y 2 ) - /3 sin(ayi) sinh(/3y 2 )) | 
= 2(/3 2 - a 2 )(B a . ) t + (a 2 + f) 2 B aS . 

□ 

In what follows, and for the sake of simplicity, we use the notation B = B a p and Ai t = (-Mq,^)^ 
if no confusion is present. 

Proposition 3.2. Let B = B a p be any mKdV breather. Then, for any fixed t € R, B satisfies 
the nonlinear stationary equation 

G[B] := B {ix) - 2(/3 2 - a 2 ){B xx + B 3 ) + {a 2 + f3 2 ) 2 B + bBB 2 x + 5B 2 B XX + ^B 5 = 0. (3.6) 
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Remark 3.1. This identity can be seen as the nonlinear, stationary equation satisfied by the 
breather profile, and therefore it is independent of time and translation parameters x±,X2 € K. 
One can compare with the soliton profile Q c (x — ct — x$), which satisfies the standard elliptic 
equation (1.51, obtained as the first variation of the H 1 Weinstein functional (1.6). 



Proof of Proposition 3.2 From (2.131 and (2.14), one has 



l.h.s. of (|3_6|) = -(B t 

= —B tx — BB X 
= -B tx + B 



2(f3 z - a z )B t + (of + (3 z yB + 5BB 2 + 5B 2 B X 



-B° 



2>B B r 



2(f3 2 - a 2 )B t + (a 2 + /3 2 ) 2 B + -B 5 



1 



-B + 2BB t - 2Mt 



-2B 2 (B t + B 3 



-B 



+ 2(/3 2 -a 2 )B t + (a 2 + /3 2 ) 2 B 

[B tx + 2M t B] + 2(fi 2 - a 2 )B t + (a 2 + f3 2 ) 2 B 



In the last line we have used (3.1). 



□ 



Corollary 3.3. Let B® - = B° g(t, x; 0, 0) be any mKdV breather as in (1.8), and x\(t), x^(t) e R 
two continuous functions, defined for all t in a given interval. Consider the modified breather 

B a . p (t,x) := B° a Jt,x;x 1 (t),x 2 (t)), (cf. ([TIT])). 

Then B a p satisfies (3.6), for all t in the considered interval. 

Proof. A direct consequence of the invariance of the equation (3.6) under spatial translations. 
Note that (3.6) is satisfied even if B a p is not an exact solution of (1.1). □ 



4. Spectral analysis 

Let z = z(x) be a function to be specified in the following lines. Let B — B a ^ be any breather 
solution, with shift parameters x\, x%. Let us introduce the following fourth order linear operator: 

C[z](x- 1) := z {ix) (x) - 2((3 2 - a 2 )z xx (x) + (a 2 + (3 2 ) 2 z(x) + 5B 2 z xx (x) + \{)BB x z x (x) 

+ [5B 2 + 10BB XX + yB 4 - 6(/3 2 - a 2 )B 2 } z(x). (4.1) 

In this section we describe the spectrum of this operator. More precisely, our main purpose is to 
find a suitable coercivity property, independently of the nature of scaling parameters. The main 
result of this section is contained in Proposition |4.11| Part of the analysis carried out in this 
section has been previously introduced by Lax [5T], and Maddocks and Sachs 22J, so we follow 
their arguments adapted to the breather case, sketching several proofs. 

Lemma 4.1. C is a linear, unbounded operator in L 2 (R), with dense domain i? 4 (]R). Moreover, 
C is self-adjoint. 

It is a surprising fact that C is actually self-adjoint, due to the non constant terms appearing in 
the definition of C. From standard spectral theory of unbounded operators with rapidly decaying 
coefficients, it is enough to prove that C* = C in H 4 (M.) 2 . 

Proof. Let z,w € i? 4 (M). Integrating by parts, one has 

/ wC[z] = [ w[z (4x) - 2(fi 2 - a 2 )z xx + (a 2 + (3 2 ) 2 z + 5B 2 z xx + lOBB x z x ] 

JR JR 

+ / [5Bl + 10BB xx + ] ^-B i -6(l3 2 -a 2 )B 2 ]zw 
Jr 2 

= f [w {ix) - 2((3 2 - a 2 )w xx + (a 2 + [3 2 ) 2 w + (5B 2 w) xx - (10BB x w) x ] z 

JR 

+ I [bBl + WBB^ + ^-B 4 - -6((3 2 - a 2 )B 2 ]zw = [ C[w]z. 
Jr 2 J R 
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Finally, it is clear that D(C*) can be identified with D(C) = H A {M). □ 

A consequence of the result above is the fact that the spectrum of C is real- valued. Furthermore, 
the following result describes the continuous spectrum of C. 

Lemma 4.2. Let a, j3 > 0. The operator £ is a compact perturbation of the constant coefficients 
operator 

C [z] := z [Ax) - 2(/3 2 - a 2 )z xx + {a 2 + p 2 ) 2 z. 

In particular, the continuous spectrum of C is the closed interval [(a 2 + /3 2 ) 2 ,+oo) in the case 
P > a, and [4a 2 /3 2 , +00) in the case /3 < a. No embedded eigenvalues are contained in this region. 

Proof. This result is a consequence of the Weyl Theorem on continuous spectrum. Let us note that 
the nonexistence of embedded eigenvalues (or resonances) is consequence of the rapidly decreasing 
character of the potentials involved in the definition of C. □ 

Remark 4.1. Note that the condition a — (3 is equivalent to the identity dpE[B] = 0. Solitons do 
not satisfy this last indentity. 

We introduce now two directions associated to spatial translations. Let B a ^p as defined in 
( fl~IT| ). We define 

Bx(t;x-L,x 2 ) := d Xl B at p(t;xi,x 2 ), and B 2 (t; x\, x 2 ) := d X2 B a ^(t; x\, x 2 ). (4.2) 

It is clear that, for all t £ K a, (3 > and Xi,x 2 £ R, both B\ and B 2 are real- valued functions 
in the Schwartz class, exponentially decreasing in space. Moreover, it is not difficult to see that 
they are linearly independent as functions of the x-variable, for all time t fixed. 

Lemma 4.3. For each t £ M, one has 

ker£ = span {Bi(t; x\, x 2 ), B 2 (t; Xi, x 2 )}. 



Proof. From Proposition 3.2 one has that d Xl G[B] = d X2 G[B] = 0. Writing down these identities, 
we obtain 

£[Ei](t; 11,3:2) = C[B 2 ]{t- x u x 2 ) = 0, (4.3) 



with C the linearized operator defined in (4.1) and B\,B 2 defined in (4.2 1. A direct analysis 



involving ordinary differential equations shows that the null space of Cq is spawned by functions 
of the type 

e ±l3x cos(ax), e ±0x sin(ax), a,/3>0, 

(note that this set is linearly independent). Among these four functions, there are only two L 2 - 
integrable ones in the semi-infinite line [0, +00). Therefore, the null space of C\jjim\ is spanned 



by at most two L 2 -functions. Finally, comparing with (4.3 1, we have the desired conclusion. □ 



We consider now the natural modes associated to the scaling parameters, which are the best 
candidates to generate negative directions for the related quadratic form defined by C. Recall 
the definitions of A a B a p and ApB a p introduced in (2.9 1. For these two directions, one has the 
following 

Lemma 4.4. Let B = B aj p be any mKdV breather. Consider the scaling directions A a B and 



ApB introduced in (2.9). Then 

J A a B C[A a B] = 32a 2 /3 > 0, (4.4) 



and 



J ApB C[ApB] = —32a 2 f3 < 0. (4.5) 
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Proof. From (3.2 1, we get after derivation with respect to a and /?, 

C[A a B] = -4a[B xx + B 3 + (a 2 + /3 2 )B], C[A p B] = 4(3[B XX + B 3 - (a 2 + (3 2 )B\. 
We deal with the first identity above. Note that from ( |2.10[ ), ( |l~3| ) and ( |2.21[ ), 

/ A a BC[A a B] = -4a [ [B xx + B 3 + {a 2 + (3 2 )B]A a B = Aad a E[B] > 0. 

JR JR 

This last identity proves (4.4). Following a similar analysis, and since E[Q] = — |JVf[Q] = — § (cf. 
(f2Tl|)-([2T2|>), one has from ( |2.11[ ) and ( |2.21[ ), 

/ A^B£[A^B] = 4/3 / + 5 3 - (a 2 + [3 2 )B]A p B 

JR JR 

= -4/3d E[B] - 16^(a 2 + p 2 ) 
= 24/3(/5 2 - a 2 )E[Q] - 16/3(a 2 + /3 2 ) = -32a 2 / 



< 0. 



Therefore, (4.5) is proved. 



□ 



A direct consequence of the identities above and Corollary |2.2| is the following result: 

Corollary 4.5. With the notation of Lemma \4-4\ let 

= aApB + /3A a B 
' 8a/3(a 2 + /3 2 ) ' 
Then Bq is Schwartz and satisfies £[Bq] = —B, 

1 



B B = 



2/3(a 2 + /? 2 ) 



> 



0, and I [ B a C[B Q ] = - 

Z JR 



1 



4/3(a 2 + /3 2 ) 



< 0. 



(4.6) 



(4.7) 



Remark 4.2. In other words, Bo is also a negative direction. Moreover, it is not orthogonal to the 
breather itself. Note additionally that the constant involved in (4.7) is independent of time. 



It turns out that the most important consequence of (4.4) is the fact that C possesses, for all 
time, only one negative eigenvalue. Indeed, in order to prove that result, we follow the Greenberg 
and Maddocks-Sachs strategy applied this time to the linear, oscillatory operator C. More 

specifically, we will use the following 

Lemma 4.6 (Uniqueness criterium, see also [H1I22])- Let B = B a ^ be any mKdV breather, and 
B\,Bi the corresponding kernel of the operator C Then C has 

^2 dim ker W[B X , B 2 ] {t; x) 

x£R 

negative eigenvalues, counting multiplicity. Here, W is the Wronskian matrix of the functions B\ 
and B2, 

B\ B 2 

{B\) x {B 2 ) x 



W[B 1 ,B 2 ](t;x) := 



(t,x). 



(4.8) 



Proof. This result is essentially contained in [TH Theorem 2.2], where the finite interval case was 
considered. As shown in several articles (see e.g. (32 US]), the extension to the real line is direct 
and does not require additional efforts. We skip the details. □ 



In what follows, we compute the Wronskian (4.8). Contrary to the 2-soliton case, where the 
decoupling of both solitons at infinity simplifies the proof, here we have carried the computations 
by hand, because of the coupled character of the breather. The surprising fact is the following 



greatly simplified expression for (4. 



Lemma 4.7. Let B = -B a ,/3 be any mKdV breather, and B\, B 2 the corresponding kernel elements 
defined in (4.2). Then 

16a 3 [3 3 (a 2 + (3 2 )[a sinh(2j8j/ 2 ) - ^sin(2a|/i)] 



detW[B 1 ,B 2 ](t;x) 



(a 2 + (3 2 +a 2 cosh(2(3y 2 ) - (3 2 cos(2ay 1 )f 



(4.9) 
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Remark 4.3. Since the computation of (4.9 1 involves only partial derivatives on the x-variable, the 
result above is still valid for the case of breathers with parameters x%,X2 depending on time. We 
skip the details. 

Proof. We start with a very useful simplification. We claim that 



detW[B u B 2 ]{t;x) = -2(a 2 + /3 2 ) / {B 2 2 (t, s) - B u {t, s)B 22 {t, s))ds 



(4.10) 



with B = B(t,x;xi,x 2 ) defined in (2.12|, and Bj = d Xj B. Let us assume this property. Using 
(2.12), we compute each term above. First of all, 



B, = 



Similarly, 



and 



Then, 



B n 
B\ 2 

B22 



B 12 



2\/2a 2 /3 cos(ayi) cosh(/3y 2 ) 
a 2 cosh 2 (/3y 2 ) + P 2 sin 2 (ay \) ' 

— 2v / 2a 3 /3sin(ayi) cosh(/3y 2 ) 
a 2 cosh 2 (/?y 2 ) + P 2 sin 2 (ayi) 

2V2a 2 /3 2 cos(ayi) sinh(/3y 2 ) 
a 2 cosh 2 (/3y 2 ) + (3 2 sin 2 (ayi ) 

— 2\/2a/3 3 sin(ayi) cosh(/3y 2 ) " 
a 2 cosh 2 (/3y 2 ) + /3 2 sin 2 (ayi) 

8a 4 /3 4 cos 2 (ayi) sinh 2 (£y 2 ) 
(a 2 cosh 2 (/3y 2 ) + /3 2 sin 2 (ayi)) 2 * 



B 2 = 



— 2v / 2a/? 2 sin(ayi) sinh(/3y 2 ) 
a 2 cosh 2 (/3y 2 ) + (3 2 sin 2 (ayi) 



2^ 2 cos 2 (ayi) 
a 2 cosh 2 (/3y 2 ) + (3 2 sin 2 (ayi) 

2a 2 cosh 2 (,Sy 2 ) 
a 2 cosh 2 (/3y 2 ) + (3 2 sin 2 (ayi) 

2a 2 sinh 2 (/3y 2 ) 



1 - 



4a 2 cosh 2 (/3y 2 ) 



a 2 cosh (f3y 2 ) + ft 2 sin (ayi) 



4a 4 cosh 4 (/3y 2 ) 



a 2 cosh (/3y 2 ) + f3 2 sin (ayi) (a 2 cosh (/3y 2 ) + (3 2 sin (ayi)) 



and 



—B\\B 22 — 



-8a 4 /3 4 sin^ayi) cosli (/3y 2 ) 
(a 2 cosh 2 (/3y 2 ) + /3 2 sin 2 (ayi)) 2 * 

2/3 2 cos 2 (ayi) - 2a 2 sinh 2 (/3y 2 ) 4a 2 /? 2 cos 2 (ayi) sinh 2 (/3y 2 ) 



a 2 cosh 2 (/3y 2 ) + /3 2 sin 2 (ayi) (a 2 cosh 2 (/3y 2 ) + /3 2 sin 2 (ayi)) 2 
Adding both terms we obtain, after some simplifications, 



^12 ~ BnB 22 — 



cos 2 (ayi) sinh 2 (/3y 2 ) - sin 2 (ayi) cosh 2 (/3y 2 ) 



(a 2 cosh 2 (/3y 2 ) + (3 2 sin 2 (ayi)) 2 

2sin 2 (ayi) cosh 2 (/3y 2 )[a 2 sinh 2 (/3y 2 ) - j3 2 cos 2 (ayi)] 



+ 



y 2 cosh (/3y 2 ) + P 2 sin (ayi) 



8a*^*fci(t,x) 



(a 2 cosh 2 (^y 2 ) + /? 2 sin 2 (ayi)) 3 ' 



with 



ki{t,x) := a sinh (/3y 2 )cosh (/3y 2 ) - a sin (ayi) cosh (,3y 2 ) 
- /3 2 sin 2 (ayi) cosh 2 (/3y 2 ) - /3 2 sin 2 (ayi) cos 2 (ayi). 



Using double angle formulas, as in ( |2.5[ )-(2.6), we get 

B12 ~ B\\B 22 



16a 4 /3 4 fc 2 
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where 



k 2 (t, x) := a 2 sinh 2 (2/3y 2 ) - P 2 sin 2 (2ayi) 



(a 2 + (3 2 ){l + cosh(2/% 2 ) - cos^a^) - cos(2ayi) cosh(2/3y 2 )). 



and g a ^ was denned in (2.17|. The last steps of the proof are the following: since 

h(t,x) := 8a 3 f3 3 [f3sm(2a yi ) - a sinh(2/?y 2 )] 

satisfies 



and 

we finally get 



{k 3 ) x = 16a 4 /? 4 [cos(2oyi) - cosh(2/% 2 )], 

S2 5 6 _ (kz{t,xy 



B\ 2 — B\\B 2 



Now, with regard to (4.10), we integrate in space, to obtain 

W[B 1 ,B 2 ] = -2(a 2 + f3 2 )'^ 1 

as desired. 



(4.91, 



We prove now (4.10). From (2.13), taking derivative with respect to x\ and x 2 , we get 
(Bi)xx + (Bi)t + 35 2 Bi = 0, (Ba)xx + (S 2 ) t + 3B 2 B 2 = 0. 



(4.11) 



Multiplying the first equation above by B 2 and the second by —B\, and adding both equations, 
we obtain 

{Bi) xx B 2 — {B 2 ) xx Bi + (Bi) t B 2 — {B 2 ) t Bi = 0, 

that is, 

{{B{) X B 2 - {B 2 ) x Bi) x = {B 2 ) t B x - (Bi) t B 2 . (4.12) 
On the other hand, since we are working with smooth functions, one has B = B\ + B 2 , 

B\ = B\\ + Bi 2 , B 2 = B\ 2 + B 22 , 

and 

(.Bi) t = SB n + jB 12 , {B 2 ) t = 6B 12 + 1 B 22 . 
Replacing in (4.12), we get 

{{B X ) X B 2 - (B 2 ) X B!) X = (5 - j)(B 2 2 - B n B 22 ). 

Since S = a 2 — 3/3 2 and 7 = 3a 2 — f3 2 , substituting above and integrating in space, we obtain the 
desired conclusion. The proof is complete. □ 



Proposition 4.8. The operator £ defined in (4.1) has a unique negative eigenvalue — Aq < 0, 
of multiplicity one. Moreover, Xq = \o(a, /3,Xi,x 2 ,t) depends continuously on its corresponding 
parameters. 



Proof. We compute the determinant (4.8) required by Lemma 4.6 From Lemma 4.7 after a 
standard translation argument, we just need to consider the behavior of the function 

/(itt) =ft, a ,p,xM :=asinh(2/3y 2 )-/3sin(2a(y 2 + (5- 7 )t + i 2 )), (4.13) 
for x 2 := xi - x 2 G R, and d — 7 = -2(a 2 + /3 2 ). 

A simple argument shows that for y 2 € R such that | sinh(2/?y 2 )| > / has no root. Moreover, 
there exists Rq = Ro(a, (3) > such that, for all y 2 > Rq one has f(y 2 ) > and for all y 2 < —Ro, 
/(2/2) < 0. Therefore, since / is continuous, there is a root yo = yoit, a, j3, x 2 ) € [— i?o,-Ro] f° r /• 
Additionally, if y 2 7^ 0, 

f'(y 2 ) = 2a/3[cosh(2/?2/ 2 ) - cos(2a(y 2 - 2(a 2 + f3 2 )t + x 2 ))} > 0, 



Miguel A. Alcjo and Claudio Munoz 



17 



by simple inspection. Therefore, if j/ 7^ then it is unique and then 

^2 dimker W[B 1 ,B 2 ](t;x) = dimker W[B U B 2 ](t; y - jt - x 2 ) = 1, 

since B\ or (B\) x are not zero at that time. Indeed, it is enough to show that W[B\ 1 B2](t, x) is 
not identically zero, then dimker B2] < 2. In order to prove this fact, note that from (4.11 1 

B\ solves, for t, Xi,X2 S K fixed, a second order linear ODE with source term —(B) t . Therefore, 
by standard well-posedness results, both B\ and [B\) x cannot be identically zero at the same 
point. 

Now, let us assume that 1/2 = is a zero of /. We give a different proof of the same result 
proved above. From (4.13), t = t^ must satisfy the condition 

-2(a 2 + p 2 )t + x 2 = ^, fceZ, 
la 

(compare with (1.13)). In terms of the variables y\ and y 2 , one has y\ = fc € Z, and y 2 = 0. 



Recall that 

B t (t,x) = -2v / 2a 2 /3 



asm(ayi) cosh(,% 2 ) + P cos(m/i) sinh(/?y 2 ) 
a 2 cosh 2 (/3y 2 ) + P 2 sin 2 (ayi) 

acos(ayi) cosh(/3y 2 ) - P sin(ayi) sinh(/3y 2 )] 



+ 2/3 sin(ayi) cos(ayi) 



Replacing in (4.14), we get 

Bi{t k ,x) 



{a 2 cosh 2 (/3y 2 ) + P 2 sin 2 (cn/i)) 2 
2y2a 3 /3sin(|7r) 



(4.14) 



a 2 + (3 2 sin 2 (|7r) 



^ 0, fc odd, 
0, fc even. 



For the first case above we conclude as in the previous one. Finally, if i& satisfies 



-2(a 2 + p 2 )t + x 2 



fcrr 



fc G Z, 



one has y\ = fcrr, y 2 = and Bi(t,x) = 0, but from (4.14), after a direct computation, (Bi) x is 
given now by the quantity 



1 



P 2 lh ,W 2 a 2 a 

— ^ cos(fc7r) H -: — COS (fcTT) 

a 4 



7^0 



for all fc <E Z. Therefore X^gk dimker W[Si, -B 2 ](ifc; x) = 1. In conclusion, for all t £ M, £ has 
just one negative eigenvalue, of multiplicity one. □ 

Corollary 4.9. There exists a continuous function f = /o(a,/3), well-defined for all a, P > 0, 
and such that 

-A 2 < -f o (a,0) <0, 

/or a/Z a, P > 0, and t,x\,x 2 S K. 

Proof. This is a consequence of the translation invariance and the fact that Ao is a continuous, 
positive function only depending on a, P and x\ := (5 — j)t + (x± — x 2 ), periodic in x\ (and then 
uniformly positive with respect to x\). □ 

Remark A A. Note that the result above is not clear if we allow a, P depending on time, as in [T5] . 
Since we do not require any kind of modulation on a and P, we can easily conclude in the previous 
result. 

Let z E H 2 (R), and B = B a .p be any mKdV breather. Let us consider the quadratic form 
associated to C: 



Q[z] 



zC[z] 



z 2 x + 2(p 2 -a 2 ) / z x + (a 2 + p 2 ) 2 / z 



2,2 



B z z 



2„2 



Biz 



10 / BB TT z' 



15 
~2 



B 4 z 2 - 6(P 2 - a 2 ) 



B 2 z 2 . 



(4.15) 
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Remark 4.5. From the definition of Q and Lemma 4.3 it is clear that <2[£>i] = Q[-B 2 ] 



0. 



Moreover, inequality (4.4) means that A a B is actually a positive direction for Q, a completely 
unexpected result. Additionally, from (4.15) Q is bounded below, namely 

Q[z] > -c Q ^||z||^ 2(R) , 

for some positive constant c a ^ depending on a and /3 only. 

Let -B-i G 5\{0} be an eigenfunction associated to the unique negative eigenvalue of the 
operator £, as stated in Proposition 4.8 We assume that £?_i has unit L 2 -n orm , so £?_i is now 
unique. In particular, one has £[-B_i] = — Aq-E>_i. It is clear from Proposition 4.8 and Lemma 
that the following result holds. 

Lemma 4.10. The eigenvalue zero is isolated. Moreover, there exists a continuous function 
vq = Voipt, /3), well-defined and positive for all a, /3 > and such that, for all zq £ H 2 (M.) satisfying 



4.3 



z B_ 1 = / z B 1 = / z B 2 = 



the 



Q[zo] > vo II z \\ 2 H 2 



(4.16) 



(4.17) 



Proof. The isolatedness of the zero eigenvalue is a direct consequence of standard elliptic estimates 
for the eigenvalue problem associated to £, corresponding uniform convergence on compact subsets 
of R, and the non degeneracy of the kernel associated to C. 



On the other hand, the existence of a positive constant vq = i'o(cv,(3,xi,X2,t) such that (4.17) 
is satisfied is now clear. Moreover, this constant is periodic in x\, continuous in all its variables, 
and satisfies, via translation invariance, the identity 

v Q (a, (3, xi, x 2 , t) = v {a, j3, xi) > 0, x x := (6 - j)t + {x\ - x 2 ), 

with vq continuous in all its variables. Thanks to the periodic character of the variable x%, we 
obtain a uniform, positive bound independent of Xi,x 2 and t, still denoted i>q. The proof is 
complete. □ 

It turns out that 2?_i is hard to manipulate; we need a more tractable version of the previous 
result. 

Proposition 4.11. Let B = B a ^ be any mKdV breather, and B\,B 2 the corresponding kernel 
of the associated operator C. There exists [1q > 0, depending on a, ft only, such that, for any 
z e H 2 (R) satisfying 



B x z = B 2 z = 0, 



has 



Q[z]>M\ 2 h* 



UfzB) 



(4.18) 



(4.19) 



Proof. This is a standard result, but we include it for the sake of completeness. Indeed, it is enough 
to prove that, under the conditions (4.18) and the additional orthogonality condition J R zB = 0, 
one has 

Q[z] > Mo||z|Ik2 (r) . 

In what follows we prove that we can replace by the breather B in Lemma 4.10| and the result 
essentially does not change. Indeed, note that from (4.6), the function Bo satisfies £[-Bo] = —B, 
and from (4.7), 



f B B = - f B a C[B a ] = -Q[B ] > 0. 
Jr jr 



(4.20) 



The next step is to decompose z and Bo in span(B_i, B\, B 2 ) and the corresponding orthogonal 
subspace. One has 



z + mB_ 1 , Bo 



nB_ 



PiBi+p 2 B 2 , m,n,p 1 ,p 2 e 
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where 

f zB-i = [ zB x = [ zB 2 = [ 60B-1 = / b B 1 = [ b B 2 = 0. 

JR il JR JR JR JR 

Note in addition that 

f B- X B X = [ B-tBa = 0. 

JR JR 

From here and the previous identities we have 

Q[z] = (Cz- mAo5_i)(2 + mB_i) = Q[z] - m 2 \ 2 . (4.21) 

jr 

Now, since £[B ] = — B, one has 

0= / zB = - zC[B ]= / C[z + mB-i]B 
Jr Jr Jr 

= I ' (C{z}-m\lB_ 1 ){b +nB_ 1 +p x B 1 +<p 2 B 2 ) = f C[z]b - mnX 2 . (4. 
Jr Jr 

On the other hand, from Corollary 14.51 



22) 



Bq£[B q 



{b + nB_!)(£[6o] - raA£B_i) = -Q[b ] + n 2 X 2 . (4.23) 



Replacing (4.22) and (4.231 into (4.21), we get 



Q[z] = Q[z] 



( f C[z]b 
v Jr 



BqB + Q[b ] 



(4.24) 



Note that from (4.20) and (4.17) both quantities in the denominator are positive. Additionally, 
note that if z = Xbo, with A 7^ 0, then 

( f C[z]b ) 2 = Q[z]Q[b }. 

K JR ' 

In particular, if z = Xbo, 

C[z}bo" 2 



I Bo 

Jr 



B + Q[b ] 



< aQ[z], < a < 1. 



(4.25) 



In the general case, using the orthogonal decomposition induced by the scalar product (£•, - ) L 2 on 
span(i3_ 1 , B 1: B 2 ), we get the same conclusion as before. Therefore, we have proved (4.25) for all 
possible z. 

Finally, replacing in ( |4.24| ) and ( |4.2l| ), Q[z] > (1 - a)Q[z] > 0, and Q[z] > m 2 A§. We have, for 
some C > 0, 

Q[z] > (1 - a)Q[z] > i(l - a)Q[z] + (1 - a)™ 2 A 2 

> ^(2PII^ (R) +2m 2 || J B_ 1 || 2 ff2(R) ) > ^||z||^ (K) . 



□ 
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5. Lyapunov functional 



In this section we introduce a new Lyapunov functional for equation (1.1), which will be well- 
defined at the natural H 2 level. 

Indeed, let uo € H 2 (R) and let u = u(t) £ H 2 (M.) be the corresponding local in time solution of 
the Cauchy problem associated to (1.1), with initial condition u(0) — uo (cf. [H]). Let us define 
the iJ 2 -functional 

F[u](t) := i J u 2 xx (t, x)dx -^-J u 2 (t, x)ul(t, x)dx + ~J u 6 (t, x)dx. (5.1) 

Lemma 5.1. Given u local H 2 -solution of (1.1) with initial data Uq, the functional F[u](t) is a 
conserved quantity. In particular, u is a global-in-time H 2 -solution. 

The existence of this last conserved quantity is a deep consequence of the integrability property. 
In particular, it is not present in a general, non-integrable gKdV equation. The verification of 
Lemma |5.1| is a direct computation. 



Using the functional F[u] (5.1 1, we introduce a new Lyapunov functional specifically related to 
the breather solution. Let B = B a ^ be a mKdV breather, and t € K, and M[u] and E[u] given 
in ((Of, (JOf. We define 

H[u](t) := F[u](t) + 2{P 2 - a 2 )E[u](t) + (a 2 + /3 2 ) 2 M[u](t), a, (3 scaling parameters. (5.2) 

It is clear that H[u] represents a real- valued conserved quantity, well-defined for i? 2 -solutions of 
(1.1). Moreover, one has the following 

Lemma 5.2. Let z € H 2 (R) be any function with sufficiently small H -norm, and B — B a p be 
any breather solution. Then, for all t€ R, one has 



H[B + z] - H[B] = ^Q[z] +Af[z] 



(5.3) 



with Q being the quadratic form defined in (4.15), and Af[z] satisfying \Af[z]\ < if||z|||p 
Proof. We compute: 



H[B + z]=\ I (B + z) 2 xx -l I 

z JR z JR 



(B + z) 2 (B + z) 2 x + - / (B + z) 



1 



+ (f3 2 - a 2 ) (B + z) 2 x --(?- a 2 ) (B + zf + -(a 2 + (3 2 ) 2 (B + z) 



1 



z JR 



Bz 



2 Tj2 



B B 



B b 



1 



B 



f [B {4x) - 2(/3 2 - a 2 )(B xx + B 3 ) + (a 2 + f3 2 ) 2 B + 5B 2 X B + 5B 2 B XX + ^B 5 ] z 

JR 1 



zi x + 2{p 2 -a 2 ) / z 2 + (a 2 +n 



?2\2 



z z - 5 / B z z 



?2 J2 



Btz 1 + 10 / BB xx z 



15 
T 



B A z 2 - 6(/3 2 - a 2 ) / B 2 z 



Bz 



.5 ' xx ^ ^ 



15 

T 



2 2 i 
Z ZZ. + 



Bz 2 x z 
Bz 5 + 



B 3 z 3 - 2(/3 2 - a' 
I ' B 2 z±-hp 2 -a 2 ) ( z'-\\ 

JR L JR Z, J R 

Therefore, we have the decomposition 

H[B + z] = H[B] + [ G[B]z(t) + ~Q[z]+Af[z], 

JR * 
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where Q is defined in (4.15), and 



G[B] = B (ix) - 2(/3 2 - a 2 )(B xx + B 3 ) + (a 2 + f3 2 ) 2 B + 5B 2 X B + 5B 2 B XX + ^B 5 . 



From Proposition 3.2 one has G[B] = 0. Finally, the term N[z] is given by 



N[z] :=5 / B 3 z 3 - 2((3 2 - a 2 ) [ Bz 3 + \ j B xx z 3 - 5 / Bz 2 x z 

JR JR 3 J R J R 

+ T / B*z±-\{p 2 -a 2 ) \z A ~\\ + U Bt + \ I z\ 

4 JR 1 JR ^ JR z JR 4 JR 

Therefore, from direct estimates one has JV[z] = 0(\\z\\^j 2 m\), as desired. □ 

The previous Lemma is the key step to the proof of the main result of this paper, that we 
develop in the next section. 



6. Proof of the Main Theorem 



In this section we prove a detailed version of Theorem |1.2| 

Theorem 6.1 (if 2 -stability of mKdV breathers). Let a,/3>0. There exist parameters r] ,A , 
such that the following holds. Consider uq € H 2 (M), and assume that there exists r\ G (0, ryo) such 
that 

Iko - ^ Q ^(0;0,0)|| ff2(R) < r?. (6.1) 
Then there exist x\{t), X2{t) S K such that the solution u(t) of the Cauchy problem for the mKdV 



equation (1.1), with initial data Uq, satisfies 



snp \\u(t) - B ai p(t;x 1 {t),x 2 {t))\\ H2{R) < A r), (6.2) 



with 

sup K(i)| + \x' 2 (t)\ <KA oV , (6.3) 

teR 

for some constant K > 0. 



Remark 6.1. The initial condition (6.1) can be replaced by any initial breather profile of the 
form B := B a p{tQ\ x\, x^), with t[j^x\,x^ £ R, thanks to the invariance of the equation under 
translations in time and space. In addition, a similar result is available for the negative breather 



-B ai p, which is also a solution of (1.1 ) 



Proof of Theorem\6J\ Let u Q e H 2 (R) satisfying (6.1 ), and let u £ C(R; H 2 (R)) be the associated 



solution of the Cauchy problem ( |1.1[ ), with initial data it(0) = uq- In what follows, we denote 

B = B aj p, 

if no confusion arises. 

We prove the theorem only for positive times, since the negative time case is completely anal- 
ogous. From the continuity of the mKdV flow for H 2 (M.) data, there exists a time T > and 
continuous parameters xi(t),x-z(t) £ R, defined for all t £ [0, To], and such that the solution u(t) 
of the Cauchy problem for the mKdV equation ( |1.1[ ), with initial data Uo, satisfies 

sup \\u{t)-B(t;x 1 (t),x 2 (t))\\ <2 V . (6.4) 
te[o,T ] 1 ' 

The idea is to prove that To = +oo. In order to do this, let K* > 2 be a constant, to be fixed 
later. Let us suppose, by contradiction, that the maximal time of stability T* , namely 

T* := sup jr > | for all t e [0, T], there exist xi(t), x 2 (t) £ R such that 

sup ||u(t)-B(t;si(t),jc 2 (t))||^ <K*r)\, (6.5) 
te[o,T] v ; ' 
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is finite. It is clear from (6.4) that T* is a well-denned quantity. Our idea is to find a suitable 
contradiction to the assumption T* < +00. 

By taking Sq smaller, if necessary, we can apply a well known theory of modulation for the 
solution u(t). 

Lemma 6.2 (Modulation). There exists rja > such that, for all rj £ (0,770), the following holds. 
There exist C 1 functions Xi(t), x 2 (t) G R, defined for all t S [0,T*], and such that 



B 2 (t;x 1 (t),x 2 (t))z(t) =0. 



z{t) :=u(t)-B(t), B(t,x) :=B aJ (t,x;x 1 {t),x 2 {t)) 
satisfies, for t £ [0,T*], 

J B 1 {t;x 1 {t),x 2 {t))z(t) - 

Moreover, one has 

\\z{t)\\ H 2 m + \x[(t)\ + \x' 2 (t)\ < KK* V , ||2(0)|U2 (K) < K V 
for some constant K > 0, independent of K* . 

Proof. The proof of this result is a classical application of the Implicit Function Theorem. Let 



(6.6) 
(6.7) 
(6.8) 



Jj(u(£),xi,a;2) 



(it(i, x) — B(t, x;xi,x 2 ))Bj(t,x;xi,x 2 )dx, j = 1,2. 



It is clear that Jj(B(t; X\, x 2 ), x\, x 2 ) = 0, for all £1,2:2 € ^. On the other hand, one has for 
d Xk Jj{u(t),xi,x 2 ) = - I B k (t,x;0,Q)B j (t,x;0 ) 0)dx. 

(B(t),0,0) 

Let J be the 2x2 matrix with components Jj^ := {d Xk Jj)j i k=i. 2 - From the identity above, one 
has 

det J = — 



B( / B}-( B 1 B 2 y (t;0,0), 
ii Jr j 

which is different from zero from the Cauchy-Schwarz inequality and the fact that B\ and B 2 are 
not parallel for all time. Therefore, in a small H 2 neighborhood of B(t; 0, 0), t € [0, T*] (given by 
the definition of (6.5|), it is possible to write the decomposition (6.6|-(|6?7|. 



Now we look at the bounds (6.8). The first bounds are consequence of the decomposition itself 
and the equations satisfied by the derivatives of the scaling parameters, after taking time derivative 
in (6.7) and using that det J 7^ 0. The last bound in (6.8) is consequence of (6.1 ). □ 



Now, we apply Lemma 5.2 to the function u(t). Since z(t) defined by (6.6) is small, we get 



from (5.3) and Corollary 3.3 



U[u]{t) = H[B](t) + -Q[z](t)+N[z](t). 



(6.9) 



Note that |JV[z](t)| < /fUzft)!!^^. On the other hand, by the translation invariance in space, 

H[B){t) = H[B](0) = constant. 



Indeed, from (1.12), we have 

B(t, x; xi(t), x 2 (t)) = B(t - t (t),x - x (t)), 

for some specific to,xo- Since H involves integration in space of polynomial functions on B,B X 
and B xx , we have 

U[B(t, ■;x 1 (t),x 2 (t))} = H[B{t - t (t), ■ - x o {t);0, 0)] = U[B(t - t (t), •; 0, 0)]. 
Finally, H[B(t - t (t), •; 0, 0)] = H[B(-, ■; 0, 0)](t - t (t)). Taking time derivative, 

d t H[B{t,-;x 1 {t),x 2 {t))] =H'[B(;-;0,0)](t-t (t)) x (l-t (t)) = 0, 



Miguel A. Alcjo and Claudio Munoz 



23 



hence H[B] is constant in time. Now we compare (6.9 1 at times t = and t <T*. From Lemma 



5.1 and (5.2) we have 



Q[z](t) < Q[z](0) + K\\z(t)f Hm) + K\\z(0)\\ 3 HHU) < K\\z(0)\\ 2 H2(m + K\\z(t)\\m 



Additionally, from ( 4.18 )-( 4.19 ) applied this time to the time-dependent function z(t), which 
satisfies (6.7), we get 

2 



<K\\z(0)\\ 2 H2m +K\\z(t)\\ 3 H2m +K 



B(t)z(t) 



< Krj* + K(K*) 3 rj 3 



K 



B(t)z(t) 



(6.10) 



Conclusion of the proof. Using the conservation of mass (1.2), we have, after expanding 
u = B + z, 



< K 



fl(0)z(0) 



K\\m\\ 2 H^ m +K\w)\\ 2 ^ m 

for each t £ [0, T*l. 



B{t)z(t) 

<K(r) + (K*) 2 v 2 ), 
Replacing this last identity in ( 6.10[ ), we get 

\\*V)\\h*( 

by taking K* large enough. This last fact contradicts the definition of T* and therefore the 
stability property (6.2) holds true. Finally, (6.3) is a consequence of (6.8). □ 



<i^ 2 (i + (ir)V)<^(inV, 
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